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CT) . ABSTRACT 
> 

fS| Planck has observed the entire sky from 30 GHz to 857 GHz. The observed foreground emission contains contributions from different phases 
, of the interstellar medium (ISM). We have separated the observed Galactic emission into the different gaseous components (atomic, molecular 
. and ionised) in each of a number of Galactocentric rings. This technique provides the necessary information to study dust properties (emissivity, 
^SJ ' temperature, etc.), as well as other emission mechanisms as a function of Galactic radius. Templates are created for various Galactocentric radii 
• using velocity information from atomic (neutral hydrogen) and molecular ('^CO) observations. The ionised template is assumed to be traced by 
' ^ free-free emission as observed by WMAP, while 408 MHz emission is used to trace the synchrotron component. Gas emission not traced by the 
^— ^ above templates, namely "dark gas", as evidenced using Planck data, is included as an additional template, the first time such a component has 
] been used in this way. These templates are then correlated with each of the Planck frequency bands, as well as with higher frequency data from 
" ^ I IRAS and DIRBE along with radio data at 1.4 GHz. The emission per column density of the gas templates allows us to create distinct spectral 
^ energy distributions (SEDs) per Galactocentric ring and in each of the gaseous tracers from 1.4GHz to 25 THz (12/jm). The resulting SEDs 
• i-H allow us to explore the contribution of various emission mechanisms to the Planck signal. Apart from the thermal dust and free-free emission, we 
have probed the Galaxy for anomalous (e.g., spinning) dust as well as synchrotron emission. We find the dust opacity in the solar neighbourhood, 
;_l ' T/Nji = 0.92 + 0.05 x 10"^' cm^ at 250 yum, with no significant variation with Galactic radius, even though the dust temperature is seen to vary from 
over 25 K to under 14 K. Furthermore, we show that anomalous dust emission is present in the atomic, molecular and dark gas phases throughout 
the Galactic disk. Anomalous emission is not clearly detected in the ionised phase, as free-free emission is seen to dominate. The derived dust 
propeties associated with the dark gas phase are derived but do not allow us to reveal the nature of this phase. For all environments, the anomalous 
emission is consistent with rotation from polycyclic aromatic hydrocarbons (PAHs) and, according to our simple model, accounts for (25 ± 5)% 
(statistical) of the total emission at 30 GHz. 

Key words. Infrared: ISM, ISM: general. Galaxy: general. Radiation mechanisms: general. Radio continuum: ISM, Submillimeter: ISM 
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1. Introduction 

In order to understand our own Galaxy, it is necessary to explore 
Galactic Plane {\b\ $ 10°), where we are able to observe the 
emission coming from a large range of distances. However the 

* Corresponding author: douglas.marshall@cesr.fr 



observed emission is the sum of a large number of line-of-sight 
components, often probing very different environments. 

Several previous studies (Bloemen et al. 19861; 
Bloe menetal 1990 ; Giard et al. 1994; Sodroski et al. 199J; 
Pala dini et al.lT2007l) have separated observed integrated emis- 
sion into a number of Galactocentric radii in order to study its 
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properties as a function of Galactic position and in different 
phases of the interstellar gas (e.g., atomic, molecular and 
ionised). The radial velocity of the gas is used to separate the 
Galactic gas emission into a number of Galactocentric rings and 
then the spectral energy distribution of each ring/gas phase is 
fitted with a physical model of dust and gas emissions. These 
methods have been used successfully in previous studies to map 
out basic properties of the ISM throughout the Galaxy. 

For instance, iGiard et aH (Il994h demonstrated that poly- 
cyclic aromatic hydrocarbons (PAHs) are a ubiqu itous com- 
ponent of the interstellar medium, and iBloemen e t al. (1990) 
showed that the dust temperature decreases with distance to the 
Galactic Centre in a way which is fully consistent with an expo- 
nential decr ease of the interstell ar radiation field (ISRF) in the 
stellar disk. Sodro ski et al ] 11997) suggested that the abundance 
of large dust grains within each gas phase exhibits a gradient that 
is equivalent, within the uncertainti e s, to th e metallicity gradient 
in the Galactic disk. iPaladini et al.l (l2007h showed that the dust 
in molecular clouds appears to be heated in a significant way by 
young massive stars still embedded in their parent clouds. 

With the advent of the Planck satellite^, it is now possible 
to perform an inversion on the emission arising from the entire 
infrared, millimetre and centimetre range, making it possible to 
determine the dust and gas properties in many different envi- 
ronments in the Milky Way. The High Frequency Instrument 
(HFI) channels (100-857 GHz) allow us to properly constrain 
the big grain temperature and emissivity, while using the Low 
Frequency Instrument (LFI) channels (30-70 GHz), it is possible 
to constrain non-thermal emission mechanisms such as free-free, 
synchrotron or anomalous dust. We have used these data along 
with other ancillary data to perform a large-scale, low-resolution 
analysis of the Milky Way ISM emission. 

Knowledge of the dust emission and physical conditions 
throughout the Milky Way provides a self-consistent context for 
other Planck studies: for exa mple, the environmental eff ect on 
the properties of cold cores (iPlanck Collaborationll20TTJd) can 
be evaluated, the modelling of individual a nomalous dust regions 
can be placed in a Galaxy-wide context jPlanck Collaboration! 
1201 Idh . or the Milky Way values can simply be compared to 
other galaxies (Planck Collaboration 2011a,b). 

In Sect. |2]we describe the data used in the study, including 
both those to be inverted and those used as templates to represent 
a particular phase of the ISM. In Sect. [3] we describe how we 
optimise the separation of the gas observation into a number of 
Galactocentric rings, as well as how we report uncertainties on 
our best solution. We present the results in Sect. [H and discuss 
potential sources of uncertainty and bias in Sect. |5] Finally we 
conclude in Sect.|6] 



2. Data 

There are two categories of data used in this study, the data to be 
inverted, and data used to create templates which will be used to 
perform the inversion. 



' Planck (http://www.esa.int/Plancki is a project of the European 
Space Agency (ESA) with instruments provided by two scientific con- 
sortia funded by ESA member states (in particular the lead countries: 
France and Italy) with contributions from NASA (USA), and telescope 
reflectors provided in a collaboration between ESA and a scientific con- 
sortium led and funded by Denmark. 



2.1. Data to be inverted 

We use data from I2yum (25 THz) to 1.4 GHz in order to study 
the dust and gas properties in a self consistent manner through- 
out the Milky Way. Here we briefly present these data in order 
of decreasing (increasing) frequency (wavelength). 

The Infrared Astronomical Satelli te (IRAS), a joint projec t 
of the US, UK, and the Netherlands dNeugebauer et al.lfl"984h . 
performed a survey of 98% of the sky at four wavelengths: 12, 
25, 60, an d 100 /vm. We use a reprocesse d version of this data 
set (IRIS, iMiville-Deschenes & Lagachd r2005) which benefits 
from a better zodiacal light subtraction and from a calibration 
and zero level compatible with DIRBE, as well as from a better 
destriping. 

The DI RBE instruinent (Diffuse Infrared Background 
Experiment. lHauser et alJl998h on board the COBE satellite im- 
aged the full sky in 10 broad photometric bands from 1 to 240 fim 
with a beam of 0.7°. The bands at 140 and 240 yum allow us to 
bridge the gap between IRAS and the highest frequencies of HFI, 
thus covering the peak of the thermal dust spectral energy distri- 
bution (SED). 

The Planck data we used for this analysis are the 
DR2 LFI and HFI maps, with the cosmic microwave 
background (CMB) removed, as describ ed in d etail in 
'Zacchei .A.. Maino. D. Baccigalupi. C. et all (l20Ilh and 
Planck HFI Core TeamI (1201 1|), respectively. It consists of a set 
of 9 frequency maps with central frequencies of 28.5, 44. 1, 70.3, 
100, 143, 217, 353, 545 and 857 GHz), smoothed to a common 
resolution of 1°FWHM, at a common HEALPix ( Gorski et al.l 
,200 5) pixelization resolution of A^side=256 (~15' pixels). For 
simplicity, in the remainder of the text we will refer to the 
first three bands as 30, 44 and 70 GHz, respectively. The CMB 
component that was removed from these maps was found by the 
Needlet Internal Li near Combination (NILC) C MB extraction 
method presented in lPlanck HFI Core TeamI ( |20I ll) . 

Finally, to provide constr aints at radio frequencies we use the 
full-sky 1.4 GHz data from (lReichlll982t iReich & Reic"hlll986l: 
Reich et al. 2001). 

All of the above data sets have been smoothed with a 
Gaussian function to a common resolution of 1°FWHM and 
gridded to a common HEALPix A^side of 256. 

2.2. Templates 
2.2.1. Atomic phase 

We aim to characterise the different emission mechanisms 
in the different phases of the interstellar medium (atomic, 
molecular and ionised). For the atomic phas e we use the 
Leide n/ Argentine/Bonn Galactic Hi Survey (iKalberla et al.l 
l200l . The survey results are expressed in antenna tempera- 
ture, Ta- In order to convert this into column density, the ef- 
fect of optical depth should be taken into account, for example 
(iBinnev & Merrifield|[T998h 

A^Hi = 1-82X 10""^ X 'r, xln^l - p-j dvcmT^, (1) 

for a given spin temperature (Ts)- R e cent s tudies in the plane of 
the Milky Way (iDickev etal.ll2003l l2009h advocate the use of 
a higher spin temperature of around 250 K. This higher value is 
justified on the grounds that along any line of sight there is a mix 
of cold and warm H i gas, and so an effective spin temperature is 
needed to correct for the opacity. This will most certainly vary 
as a function of Galactic longitude. However, constraining this 
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Fig. 1. Templates used in the inversion. The top five images are the H i templates, the following five are the molecular templates, all 
of which are presented in increasing Galactocentric radius (the radii defining each ring are described in Sect. 13.2b . The final threg 
templates are the free-free, synchrotron, and dark gas templates, respectively. The units are all expressed in 10^° atoms, except for 
synchrotron, expressed in MJy sr ' . The centre and anti-centre regions have been masked as no distance information can be obtained 
from the radial velocity of the gas. 
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Fig. 2. Longitude profile of H i, CO and dark gas. The black solid 
line shows the total contribution from H i & CO. The dark gas 
represents a significant fraction of the gas column density, and 
dominates the CO outside of the molecular ring. 

variation is beyond the scope of this study. We have therefore 
assumed a constant hydrogen spin temperature of 250 K for the 
entire Galactic Plane. Varying this value in the range 150-400 K 
has only a small effect on our results and they remain consistant 
within the quoted uncertainties. 

2.2.2. Molecular phase 

We assume that the bulk of the molecular gas mass is at first or- 
der well traced by CO emission (Tielens 2005) . The CO data w e 
use come from the Composite CO Survey of lDame et al.l(l200lh . 
The column density of molecular hydrogen can be expressed as 
a function of antenna temperature (Tco) 

Nn,^2x Xco X f Tco dv. (2) 
Ji'i 

where Xco - A^(H2)/VVco is the ratio of molecular hydrogen 
column density to the velocity-integrated intensity of the '^CO 
line. We use a value for Xco which is comp atible with the re- 
cent values from the F ERMI collaboration dAbdo et al.ll2010t 
lAckermann et al.ll201 ih . Xco = 1-8 x 10^", which results from a 
study on diffuse gamma-ray emission in the Galactic plane. 

Both the atomic and molecular data contain line of sight ve- 
locity information. It is therefore possible to use this informa- 
tion to separate the emission into a set of Galactocentric rings, 
as a given line of sight velocity measurement uniquely identifies 
a Galactocentric distance, assuming circular orbits. This should 
not be confused with heliocentric distances for which there is 
a distance ambiguity corresponding to the near or far side of 
the derived Galactocentric orbit. The velocity information alone 
does not allow one to determine if the observed gas is on the 
near side of the determined orbit or on the far side. We sepa- 
rate the observed H i and CO observations each into a set of 5 
Galactocentric rings. The reasons for choosing this decomposi- 
tion is explained in Sect. l3.2l 

2.2.3. Dark gas 

It has been shown (e.g., iReach et al]ll994t iGrenier et ani2005h 
that the CO and H i do not account for all the atomic and molecu- 
lar gas in interstellar clouds: there is in addition a dark gas phase. 
This ph ase is also extensively studied in Planck Collaboration! 
(l20Tl^ . where the total gas column density, as traced by H i and 
CO observations, is correlated with the optical depth of the dust 



calculated using IRAS 100 yum along with Planck HFI data at 350 
and 500 yum (857 and 545 GHz). The correlation is linear at low 
and high column densities, but it departs from linearity at col- 
umn densities of 7.4 x 10^° < A^h < 5.0 x 10^'. This excess is 
attributed to thermal emission by dust associated with a dark gas 
phase, undetected by H i and CO measurements. 

The dark gas column density used here is defined in 
iPlanck Collaboration (.201 Ic.) as 

(\ lef 

where tq is the thermal dust optical depth, tm = 
i+axcoWcol) is the modelled dust opacity and 
(tdINwY^^ is the reference dust emissivity measured in low A^h 
regions. We use their dark gas estimate at 857 GHz as an addi- 
tional template. 

The dark gas phase represents a significant fraction of the 
total gas column density, as shown in Fig.|2] The Hi dominates 
over almost the entire longitude range, and the CO is strongest 
in the molecular ring 40° > I > -40°. Outside of the molec- 
ular ring the dark gas is stronger than the CO, except towards 
the anti-centre direction where the dark gas and CO are roughly 
comparable. The proportion of dark gas to total hydrogen col- 
umn density ranges from roughly 10% in the anti-centre direc- 
tion up to nearly 60% towards the Galactic centre. 

2.2.4. Synchrotron emission 

An additional emission mechanism which is important in the fre- 
quency range we are studying and is not traced by usual gas 
tracers is synchrotron emission, which arises from the acceler- 
ation of cosmic -ray electrons in magnetic fields. This emission 
follows approximately a power law but the spectral index most 
certainly varies on the sky and as a function of frequency. More 
than 90% of the o bserved synchrotron emission arises from a dif- 
fuse component (iBennett et al.ll2003l) and in the Galactic plane 
recent studies have shown that there is no hardening of the syn- 
chrotron spectrum tFinkbeiner 2004; Boughn & Pober 200'^, 
but that emission from spinning dust is likely to dominate be- 
tween 8 and 60 GHz. For our synchrotron templa te we hav e 
chosen the 408 MHz all-sky survey from Haslam et al.l ([T98l . 
which traces the diffuse, soft synchrotron component. 

2.2.5. Ionised phase 

To represent the ionised component we have used the Maximum 
Entropy Metho d (MEM) free-free template from the WMAP 
seven year data (iJarosik et al.ll20I ll) . The intensity map (7^) has 
be en converted to ionised hydrogen column density as described 
bv lSodroskiet al.l(ll997h : 

Nn„ = 1.2xl0'5rf%°in;^/, (4) 

for effecti ve electron den s ity n^s and electron temperature T^. 
Following ISodroski et all (figsl . we use Heff = 10 cm and 
Tf. = 8000 K. In the future, the availability of radio recombina- 
tion line (RRL) surveys will provide a more direct measurement 
of the free-free component in the Galaxy and its distribution sep- 
arated into distinct Galactocentric bins using the radial velocity 
information included in such surveys. 

As has been done for the data, the templates have all been 
smoothed to 1° FWHM and gridded to an A^side of 256. These 
templates are all displayed in Fig.[T] 
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Fig. 3. Results of the inversion for the 857 GHz band. From top to bottom are: observed emission; modelled emission; atomic 
contribution; molecular contribution; ionised contribution; dark gas contribution; synchrotron contribution; residual; and relative 
residual. All images except for the residual maps are at the same intensity scale. The centre and anti-centre regions have been 
masked, as the kinematic distance method is inapplicable to these regions. Areas where no CO observations are available have also 
been masked. 
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3. The Galactic inversion 

Component separation techniques have been widely tested in the 
framework of multi-frequency observations of the CMB, espe- 



cially in the context of Planck (see lLeach et al.ll2008l , for a re- 
view of the methods developed within the Planck collaboration). 
These methods exploit spectral and spatial correlations and inde- 
pendence in observations made at separate frequencies, as well 
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as external constraints and physical modelling, in order to dis- 
tinguish between different physical sources of emission. Most 
of them are designed to efficiently separate CMB from fore- 
grounds and possibly foregrounds into "components": (Galactic 



synchrotron, free-free and dust emissions; extra-galactic and far- 
IR point sources; Sunyaev-Zeldovich eff'ect, etc. 

The approach presented here is different from these methods 
as it tries to quantify the separate components responsible for 
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diffuse Galactic emission. Nevertheless it is not strictly speaking 
a component separation method, in the sense presented above, 
as it fits one map of the sky with spatial templates at different 
Galactocentric distances and is not exploiting the spectral signa- 
ture of each component to perform the separation. 

3. 1 . Inversion method 

We assume that each of the frequency maps (ly) can be expressed 
as a linear combination of the spatial templates described in 
Sect. l2.2l and presented in Fig.[T] Then, for each frequency map, 
we are solving for the mean emissivity of each template (in units 
of MJy sr"' / 10^" atoms) which best reproduces the observed 
data. Note that we are not taking into account any spatial varia- 
tions of the emissivity within a given template, but are assuming 
a constant value for each of our templates. The modelled emis- 
sion can be written, for any frequency map Iy(x,y), as: 

i:(x,y) = 2ti(4,(v)<(^,3') + 4M^H,(^'3')) 
+ eH„{v)NH„ix,y) + esiv)N,{x,y) 
+ ediv)Ndix,y) + Cy, (5) 

where the sum is over a number (n) of Galactocentric rings (/) 
for the atomic and molecular phases, denoted A^Hi and Nhi re- 
spectively. A^H„ , Ns and A^d denote the ionised, synchrotron and 
dark gas components, respectively, and Cy is a constant. Note 
that the last three templates do not have any rings, as we have no 
way of separating the information into different Galactocentric 
distances and so they are a function of frequency only. As such, 
the quoted e are a function of v only. 

In order to find the values for the different emissivities, e, we 
minimise the chi-squared value defined by 

where o-y is the noise in the map ly, estimated as the median 
value of the intensity for \b\ > 70°. 

In our case, the difference between model and observations 
is not at all dominated by the low noise present in the Planck 
bands, but by the simplistic model that we are adopting for the 
inversion. As such we suppose that the noise does not vary spa- 
tially. 

The solution which minimises Eq.|6]is found via matrix in- 
version. By writing the problem as 

bi = fly ■ Aij, (7) 

where b is the observed intensity map, A is an x M matrix of 
the different templates and a is the vector of the parameters. The 
solution is simply (Press et al. 1992) 

a^iA^Ay^ xA^b . (8) 

The optimal mathematical solution which solves Eqs. |5}|6] 
may not be physical; in particular, it may include negative emis- 
sivities. However, we do not have any cases where the emissivity 
is more than one sigma below zero, so remaining consistent with 
zero. As such we do not place any constraints on the parameters 
of the inversion. 

3.2. Galactocentric ring optimisation 

The observed radial velocities of H i and CO are used to sepa- 
rate the observed emission into a number of Galactocentric rings. 



assu ming circular o rbits for the gas and using the rotation curve 
from lClemensI (119851) . The choice of the distance intervals for the 
rings is not arbitrary, but is subject to a number of constraints. 
Choosing too many would create a set of rings that would be 
highly correlated and result in a degenerate solution. An opti- 
mised set of rings has been defined by minimising the correlation 
between rings, and thus minimising the degeneracy between the 
various components of the inversion. The total number of rings 
and their thickness are also driven by the requirement of a min- 
imal number of data falling into each ring in order to constrain 
it and separate its contribution from the others. Further, we have 
chosen to fix the solar circle as one of the rings, including all 
velocity channels that contain high latitude, low velocity clouds. 

The ring optimisation process consists first in dividing the 
H I and CO templates into a large number of small rings, and 
then building the covariance matrix {A^A, see Eq.[8]) associated 
with this highly discrete basis of annuli. An iterative process 
is then performed: neighbouring rings with the highest corre- 
lation are merged and the new covariance matrix is updated with 
a reduced number of annuli. The process continues, merging the 
rings which have strongest correlation between them, until a rea- 
sonable number of rings is obtained. This leads to the following 
decomposition: Rqal = [1.5, 4.2, 5.4,7.1, 9.8, 25.0]kpc, where 
the values define the ring edges. 

3.3. Estimating uncertainty 

The result of the inversion is an average gas emissivity at each 
frequency, in each gas phase and in each ring (when available). 
However, within a given gas phase the dust and gas properties 
may vary significantly and we wish to report this uncertainty. In 
other words, the residual differences between the model and the 
data are not due to noise in the data, meaning that the uncer- 
tainty on the resulting emissivities cannot be computed analyt- 
ically (i.e., the residuals are not well-modelled by the uniform 
Gaussian noise implicit in the solution, Eq.[8]l. The residuals are 
rather due to unknown differences between the templates and the 
data. An alternate method is then needed to estimate the uncer- 
tainties on the results. 

The dust properties in the northern Galaxy may differ from 
those in the southern part, and of course the variation may be 
even more localised than that. To quantify this uncertainty we 
perform jackknife tests by masking approximately 30% of the 
Galaxy in 30 random blocks, each 3°wide, and performing the 
inversion on the unmasked zones. This is repeated 500 times, 
and the rms value of the jackknife tests is used to provide the 
total variation of the emissivity within each ring and gas phase. 
We therefore provide robust estimates of the dust emission and 
provide their variation for each ring and gas phase. 

It should be stressed that this estimation does not take into 
account the errors on either the H i opacity correction nor the 
uncertainty on the assumed Xqo value for each phase globally. 

4. Properties of the ISM in the Galactic plane 

The above inversion method has been applied to 16 different fre- 
quency bands from 12 /vm (25 THz) to 1.4 GHz. The results are 
presented in table form in the appendix (Table lATI ). 

Figures[3]|4]and|5]compare our best-fit model to the data for 
three frequencies (two Planck frequencies, 857 and 30 GHz, as 
well as the 1 .4 GHz radio frequency) and the contributions of the 
different gas components to the model. The model reproduces 
the data rather well, especially when one notices the scale of the 
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residual. The contours in the relative residual are drawn every 
25%, so nearly the full zone used in the inversion has residuals 
below this level. 

The main contributors to the 857 GHz map (see Fig. [3]) are 
the atomic and molecular phases (including the dark gas) with 
a modest contribution from the ionised component. At 30 GHz 
(Fig. ID the emission is a combination of the molecular, ionised, 
dark gas and synchrotron, with very little hydrogen. As for the 
solution at 1 .4 GHz (Fig. |5]) it is the synchrotron component that 
dominates the solution, but with a significant contribution from 
the ionised component. 

4. 1 . Dust in thermal equilibrium 

We have constructed SEDs for each of our templates across all 
frequencies that we have analysed. Where the emissivity is com- 
patible with zero, we plot the 2cr upper limit as a downward 
arrow. 

We have fit the SEDs using a sum of thermal dust (big 
grains), spinning PAH molecules, free-free and synchrotron 
emission. The large dust grains in thermal equilibrium are well 
described by a modified blackbody of the form 



l)vo ' 



(9) 



where T is the dust temperature and p is the emissivity spec- 
tral index. When fitting using both T and /? as free parame- 
ters, we find a mean value of y6 - 1.8 ± 0.1 throughout the 
Galaxy. However, to i ncrease the reliability of the fit, we follow 
iPlanck Collaboration! (f2011g.h.c) and fix j6 = 1.8. We use all 
data points down to 143 GHz in the fit, but exclude 217 GHz in 
the molecular and dark gas phases due to contamination from 
CO line emission. All points are colour corrected assumin g 
v/y = constant as described in IPlanck HFI Core TeamI (1201 Ih . 
The SEDs for the atomic, molecular and ionised phases are 
shown in Figs. |6] [T] and [8l respectively. We are able to con- 
strain the dust temperature for each component and in the var- 
ious Galactocentric radii. The resulting temperature profile is 
shown in Fig |9] The temperature in the atomic phase is seen 
to decUne steadily with increasing Galactocentric distance. This 
decrease in temper ature follows the trend of the ambient ISRF 
(iMathis et al.lll983 ). shown by the solid line in Fig|9] indicating 
that the dust in the atomic phase is predominantly heated by it. 

The dust in the molecular component, however, is seen to be 
more steady, peaking in the molecular ring (3 ^ ~ 6 kpc) 
and once again outside the solar circle. This phase seems to be 
heated by the presence of star formation. The ionised phase (dot- 
ted line) is seen to be dominated by grains that are on average 
much warmer than either of the other two phase s. These results 
are i n agreement with previous Galactic analyses ('S odroski et al.l 
1199 7: Paladini et al. 2007). The temperature within the solar cir- 
cle {T - 17.6 ± O.IK) is completely in agree ment with the 
average value at high Galactic latitude (.Planck Collaboration! 
1201 lgl) . 

The dust opacity at a given wavelength is defined as t^/Nh 
where - £^/B(v, T) is the dust optical depth. The dust opac- 
ity in the atomic phase for the solar circle is equal to 0.92 ± 
0.05 X 10"2-**cm2 at 250 yum, slightly below, but in agreem ent 
with, the same value found by IPlanck Collaboration! (1201 Igh at 
high Galactic latitude. The variation of the opacity as a function 
of Galactic radius is shown in Fig.[TO] We do not detect any sig- 
nificant variation of the dust opacity as a function of Galactic 
radius, in either the atomic or molecular phases. Note, however. 
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Fig. 6. SEDs for the atomic component, each sorted in increas- 
ing Galactocentric radius. The colour of the symbols refer to the 
mission the data came from: IRAS (blue), DIRBE (yellow), HFI 
(red), LFI (green) and 1 .4 GHz (black). The SEDs also show var- 
ious fitted laws: the dotted line is the thermal dust SED, the short 
dashed line spinning dust in an atomic medium, and the solid 
line the sum of all contributions. The spinning dust contribu- 
tion is not a fit, but simply the result of using typical values in 
the model (see Sect. 14.2.31 for details). The contribution of the 
atomic component in the outer Galaxy to the 1 .4 GHz signal is 
lower than the noise level derived for the 1 .4 GHz map and so 
should not be viewed as significant. 
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Fig. 7. SEDs for the CO component, each sorted in increasing 
Galactocentric radius. The different Hnes and colours are as de- 
scribed in Fig. |6] except for the long dashed line which repre- 
sents spinning dust in a molecular medium. 



that the opacity of the molecular component is completely de- 
generate with the chosen value for Xqo- By increasing Xco, we 
decrease the opacity. Further, the appropriate spin temperature 
for the H i opacity correction may differ from the constant 250 K 
we have assumed. 

Finally, after fitting the thermal dust in the molecular 
phase, there is an obvious excess at 100 and 217 GHz due 
to the presence of CO line emission in these Planck bands 
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Fig. 8. The resulting SEDs for the free-free component. The dif- 
ferent lines and colours are as described in Fig.|6] except for the 
dash-dot line and the dash triple dot line which represent free- 
free emission and spinning dust in an ionised medium, respec- 
tively. 
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Fig. 9. Average dust temperature as a function of Galactocentric 
radius for the atomic (blue) and molecular (red) components. 
The dust temperature in the ionised phase and the dark gas phase 
are shown by the dashed line and dot-dashed line, respectively. 
The solid line shows the behaviour of the ISRF with Galactic 
radius. Horizontal bars show the range in distance over which 
the mean temperature applies. The temperature of dust mixed 
with atomic Hi decreases with radius, consistent with the fact 
that this component is heated by the ambient ISRF. Dust in the 
molecular component is heated by embedded star formation. 



dPlanck HFI Core TeamI 1201 lb . This is a clear indication that 
the inversion is working as it should, as this excess is present 
in the molecular phase. The excess with respect to the fit is 
16.2+1.0 an d 20.2+4.5 a<Kcmb/Krj km s"' at 100 and 217 GHz, 
respectively. iPlanck HFI Core TeamI ( |201 ll) find 14.2 ± 1.0 and 
44.2 ± 6.0 at 100 and 217GHz, respectively. At lOOGHz our 
value is slightly higher, which is to be expected as it is simply an 
excess with respect to '^CO and does not take into account the 
other CO isotopes present along the line of sight. The difference 
in the values at 217 GHz is not as easy to understand, but may 
come from our simple assumption of a fixed yS. 



4.2. Analysis of non-thermal SED components 

The SEDs for all components at frequencies lower than 100 GHz 
show non-thermal dust emission features. The emission mecha- 
nisms at work in this spectral range are thermal bremsstrahlung 
(free-free), synchrotron emission, as well as a third commonly 
referred to as anomalous dust emission. 
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Fig. 10. Dust opacity at 250 /iin as a function of Galactocentric 
radii for Hi (blue, diamonds) and CO (red, triangles). The 
dashed line and the dash dot line show the dust opacity for the 
ionised phase and the dark gas phase, respectively. 
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Fig. 11. The resulting SEDs for synchrotron component. The dif- 
ferent lines and colours are as described in Fig. |6l The point at 
408 MHz is not part of the inversion, but has been plotted here 
to show consistency. 



4.2.1 . Synchrotron 

The SED for the synchrotron component resulting from our in- 
version is shown in Fig.[TT] A power-law has been recovered in 
the two lowest frequency LFI points and the 1 .4 GHz radio point. 

For electrons with a power-law distribution of energies 
A'(£') oc E^P the frequency dependence of the emission is char- 
acterised by antenna temperature T{v) oc v'^ with spectral index 

jS = -{p + 3)/2, with typically fi 3. As we are not working 

in antenna temperature, but flux density, the power law takes the 
form ly oc v" with a - p -2. 

We have performed a fit of the rising synchrotron spectrum 
with a power law and obtained a spectral index a = -1.0 ± 0.1. 
This value is slightly flatter than the averag e value found be- 
tween 408 MHz a nd 23 GHz by WMAP (Du nklev etalll2009t 
iDavies et~d! I I2OO6I) . However, this value f alls within the range 
of values found in the Galaxy jBennett et a l. 2003, -0.5 < a < 
-1.1), and most likely represents the stronger power injection 
in the Galactic plane. There is no evidence in the data for any 
hardening of the synchrotron spectrum below 70 GHz. 

4.2.2. Free-free 

Free-free emission dominates the ionised hydrogen template 
(Fig [8| at frequencies below 143 GHz. This emission follows 
a power-law with a spectral index of a - -0.1 for 70GHz< 
V < 1 .4 GHz. As we have used the WMAP free-free template to 



trace the ionised medium, it is interesting to compare the total 
free-free emission that we attribute to this phase with the WMAP 
estimate. Using the power-law fit shown in Fig. |8] we can ex- 
trapolate to the frequency of our choice and, using our H-t- tem- 
plate, produce an estimate of the free-free emission in any of the 
WMAP bands. 

Performing this at 23 GHz, we recover ~ 80% of the emis- 
sion that the WMAP team find in their MEM free-free component 
separation. Using radio recombination fin es (RRLs) from the H i 
Parkes Zone of Avoidance (ZOA) survev. ^lves et al] (1201 Ol) re- 
cently showed that diffuse ionised emission can be recovered in 
the Galactic plane using RRL surveys. They then compared the 
level of free-free derived with the WMAP MEM estimate. As in 
this work, they showed it to be too high in the Galactic plane by 
at least 10%, albeit in a much smaller (8° x 8°) region. 

Looking at the molecular SEDs (Fig. |7]l, and the dark gas 
SED (Fig. [OJ, the LFI points depart from the thermal dust and 
seem to describe a free-free spectrum. It is not possible, however, 
to describe the three LFI points and the radio point using a free- 
free power law, except for the first and last molecular templates. 
Even so, all the radio points in the molecular templates are upper 
limits and the predicted free-free spectrum would be above the 
4-cr level for three of the templates, showing that the data are 
not compatible with this model. In addition, if this were indeed 
free-free emission, our estimate of the total Galactic free-free 
emission would surpass that of the WMAP MEM fr ee-free esti- 
mate b y over 30%. This seems unlikely in light of the lAlves et al.l 
(1201 oh finding discussed above. In that case, the LFI-radio points 
in the atomic, molecular and the dark gas phases are ano ma- 
lous microwave emission. Indeed, iBoughn & Pobeii (l2007h re- 
cently showed that the dominant emission mechanism at 20 GHz 
is spinning dust, and so is identified as anomalous microwave 
emission. 



4.2.3. Anomalous microwave emission 

Anomalous microwave emission, ari sing between 10 and 
100 GHz and first detected in the 1990s (iKogut et al.lll996h . can- 
not be explained using classical emission mechanisms known 
in this frequency range CBandav et al., .2003). Indeed it is too 
bright to be free-free emission, and it is not polarised as ther- 
mal dust or synchrotron emissions would be. However, it is 
correlated with dust IR emission and especially with the inter- 
stel lar PAH and very-small-grain (VSG) emission in the mid - 
IR (iCasassus et al.ll2006t lYsard et al.ll2010l IScaife et aI.ll20Toh 
Anomalous emission was first proposed by Draine & LazarianI 
(11998*) to come from rapidly rotating PAHs ("Rouan et al.l| l992h . 
the spinning dust grains. Since, it has been observed in vari- 



ous interstellar environments: dark clouds dWatson et al] 20051: 

aLlI; 



ICasassus et aL II2OO6I: ^aife et al. "20091 'Dickinson et al.ll2oTol) 



H n regio ns (iDickinson et al. 2007; Scaife et al. 2008); planetary 
nebulae (Casassus et al. 2007) and diffuse interstellar medium 
(Miville-Deschenes et al. 2008). Several processes can excite 
or damp the grains' rotation: photon emission (IR and radio); 
gas/grain interactions; formation of H2 molecules on the grain 
surface; or photoelectric emission. 

Here we investigate the anomalous component in CO- 
correlated emission. All the rings exhibit a flattening of the SED 
in the three LFI bands (Fig. |7]l, as well as upper limits on the 
emission at 1.4 GHz. This clear excess in the microwave range 
cannot be free-free emission due to the reasons discussed in 
Sect. 14.2.21 The peak frequency of this excess lies approximately 
from 20 to 60 GHz. 
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Fig. 12. Depiction of the spinning dust emission at 30 GHz in the Galactic plane in MJy sr ' (top) and as a fraction of the total30 
GHz signal (bottom). In making this image, it is assumed that the 30 GHz emission associated with the atomic, molecular and dark 
gas phases are solely due to spinning dust. Some features may be simply due to residuals in our inversion; however, this illustrative 
view highlights regions where spinning dust may be particularly strong compared to other emission mechanisms. 



We attempt to fit the e xcess with the spinning dust model 
described in Silsb ee et al.l ([2011). Models show that spinning 
dust emission is sensitive to the grain size distribution, to 
the intensity of the ISRF, and to th e density of the medium 
in which the grains are embedded (lAli-Haimoud et al.) l2009t 
lYsard & Verstraetd l2010h . The smallest grains, namely the 
PAHs, are responsible for most of the spinning dust emissiv- 
ity. Fitting the thermal dust emission with a modified blackbody, 
we estimate the rad i ation field to be clo se to the standard ISRF 
dMathis etal.1 119831; iPorteret al.ll2008l) . To describe the small 
dust grains, we use two grain populations: PAHs and VSGs. We 
adopt a size distribu tion tha t is th^siim of two log-normal com- 
ponents ( Wei ngartner & Drainell200lh . The first one is peaked 
at 3 nm (VSG), and the sec ond one at 0.6 nm (PAHs), in line 
with lCompiegne et al.l(l201 ll) . with the width of both log-normal 
components set to 0.4. Most of the spinning dust emission comes 
from the smallest grains, the second log-normal component, for 
which we assume an abundance equal to the value measured in 
the solar neighbourhood. Finally we describe the grains' envi- 
ronment as molecular gas with typical parameters («// = 350 
cm""* and Tgas = 20 K). 

The fit to the data with this simple model is surprisingly 
good, and shows that not only is anomalous emission present 
in the molecular phase throughout the Milky Way, but that spin- 
ning dust provides a compelling model with which to describe 
it. The predominance of spinning dust in the molecular phase is 
in agreement with more detailed modelling of individual anoma- 
lous emission regions (Planck Collaboration 201 Id). 

The ISRF could in theory also have an effect on this emis- 
sion mechanism, and the big grain temperature gives us an idea 
of its fluctuations: x - {T /ll .SK'f^'^, where x is dimensionless 
and equal to one in the solar neighbourhood. However this range 
of IS RF fluctuations has little effect on the spinning dus t spec- 
trum (1X5 -Hai'moud et al .l2009HYsard & Verstraetel20Toh . so we 
assume a solar value for all SEDs. The average density of the 
molecular gas may vary with radius, which would have enough 
impact to significantly alter the fit. However, we have not pro- 
ceeded with a detailed fit of this component, as there are many 
parameters and the solution would be degenerate (e.g. grain 
size distribution, gas density, ISRF intensity, as well as radia- 
tive transfer considerations). We have kept the analysis simple, 
which allows us to say that we have definitely detected anoma- 
lous microwave emission, that it is associated with the molecular 



phase, and that spinning dust provides a good model with which 
to understand the spectrum of this excess emission. Some of this 
emission could be explained by free-free emission, but it cannot 
be the dominant emission mechanism in this phase. 

To study the possibility that spinning dust emission is present 
in either the atomic or ionised phase, we have also estimated its 
contribution to these. The different environments will have an 
impact on the spinning dust emission, so we have adapted the gas 
parameters i n the model ling. For the atomic phase we have used 
the results of lHeiles & T roland (2003) and modelled the atomic 
ISM as composed of 40% cold neutral medium (CNM) and 60% 
warm neutral medium (WNM). The parameters for these phases 
are na = 30 cm"^ and Jgas = lOOK for die CNM, and ria = 
0.4 cm"^ and Tgas = 6000 K for the WNM. As can be seen in 
Fig. |6] the SEDs show that spinning dust emission is able to 
reproduce the departure from the thermal dust spectrum for v < 
70 GHz. The 1 .4 GHz point in the last atomic template is almost 
certainly an artifact. The atomic and synchrotron templates both 
describe diffuse smooth emission in the Galaxy, leading to some 
"crosstalk" between the two components. However it should be 
noted that when the emissivity found for the atomic phase in the 
outer Galaxy is multiplied by the column density of hydrogen 
in the template, the total emission is very low, lower than the 
noise level derived for the 1 .4 GHz map. As such, it should not 
be viewed as significant. 

For the ionised phase, we use the same gas parameters as in 
Sect.O namely «h = 10 cm"^ and Tg^s = 8000 K. The SED in 
Fig.[8]shows that spinning dust may also be present in this phase, 
but that it is completely dominated by free-free emission. In fact, 
the situation may be even worse, as we have assumed so lar value 
PAH abundance, whereas recently iDobler et al.l ( l2009l) showed 
that the PAH abundance in this type of region is about a factor 
of three lower. 

Assuming that the atomic, molecular and dark gas contribu- 
tions at 30 GHz are due solely to spinning dust, and that the 
ionised phase has a negligible contribution, we have produced a 
Galactic -plane image of spinning dust emi ssion. This provides a 
comple mentary image to that presented in Planck Collaborati"onl 
as it has been obtained using a different approach. The 
intensity map is shown at the top of Fig. [T2j at the bottom we 
show the fraction of the total signal at 30 GHz due to spinning 
dust. Several bright regions show up in this map, highlighting 
areas where spinning dust emission is the dominant source of 
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Fig. 13. The resulting SED for the dark gas component. The dif- 
ferent lines and colours are as described in Fig. |6] The dashed 
and short-dashed lines shows the spinning dust contribution 
from a molecular phase, and from an atomic phase, respectively. 



emission at these frequencies. Overall, spinning dust accounts 
for 25 + 5%(statistical) of the total emission at 30 GHz. This 
error estimate does not take into account any systematic uncer- 
tainty due to our simple model, which is difficult to quantify, and 
so will likely be higher. 

4.3. Dark gas 

Dark gas has been detecte d in numerous studies ( iReach et alj 
Il994t brenieretaLl l2005h as an infrared, and a correlated 
infrared/gamma-ray, excess over and above what is expected 
given the gas column density derived from H i and CO data. The 
inversion performed without this component results in positive 
and negative residuals, as the other components compete to ac- 
count for this missing gas phase. Not only does adding this com- 
ponent lower our reduced fits across all bands, but it also 
allows us to examine the nature of this dark component. 

Looking at Fi dT3l the resulting SED for this component 
looks quite similar to both the atomic and molecular SEDs. Its 
temperature does not differentiate it immediately from either of 
the two phases, and it shows a low-frequency tail beyond the 
thermal dust SED. These points look similar to the anomalous 
dust in the molecular phase (dashed line); however, the spinning 
dust model lies above the LFl points. Assuming that the spinning 
dust is in an atomic phase (dotted line), it can be seen that it does 
not predict the observed emission. It thus seems to be more of a 
mix of the two phases, with a larger influence from the molecu- 
lar phase. The different parameters discussed in Sect. 14.2. 3] could 
also be altered, but we do not have enough data to constrain such 
a fit. 

The interpretation of these results is somewhat complicated 
for several reasons. Firstly, the amount of dark gas is not ex- 
pected to be independent of Galactic position, but using a 2D 
template we are forced to assume just that. Secondly, the esti- 
mate of this gas phase was obtained at high galactic latitude and 
we are using an extrapolation into the Galactic plane. It would 
be preferable to use an independent template, derived from 
gamma-ray observations for example, in future inversion analy- 
ses. Reconstructing the line-of-sight distribution of the dark gas 
would enable a more exact determination of the properties of 
the dark gas, by co upling th e analysis with a 3D extinction map 
(e.g. lDrimmeletal . 2003; MarshaU et al.l2006l;ISale et al.l2009l) . 
This type of analysis will be considered in a future study of the 
Galactic inversion. 

In summary, the dark gas SED is consistent with molecular 
gas, most likely diffuse H2 and optically thick CO. The best fit 



at longer wavelengths comes from the spinning dust model, in a 
molecular medium. 

5. Possible sources of bias 

Our inversion algorithm determines the flux density per col- 
umn density of each template, and hence the result depends 
on the different templates themselves. Care should therefore be 
taken to use the most physically representative input for each 
phase of the ISM that we want to include in the inversion. 
Nonetheless, the jackknife tests described in Sect. l3.3l show that 
the method is self-consistent. Moreover, many different combi- 
nations for the H iH-CO Galactocentric rings has been attempted, 
and the main results presented here remain coherent with those. 
However, the results are notably sensitive to one particular in- 
put, the synchrotron template. Without it, the H i templates trace 
the Galactic synchrotron emission, as the H i observations trace 
the large-scale smooth structure of the Milky Way, similar to the 
diffuse soft synchrotron emission. This artificial component in 
the Hi templates then masks any contribution for the anoma- 
lous microwave emission. We have attempted to trace the syn- 
chrotron emission with two different synch rotron template s: the 
WMAP MCMC synchrotron template (Gold et al.' 2009); and 
the 408 MHz fuU sky survey (Haslam et al. 1982). We have de- 
cided to use the latter as it more successfully removes the syn- 
chrotron contribution from the Hi templates and allows us to 
detect anomalous emission in the atomic phase. 

5.1. Correlation between parameters and uniqueness of the 
solution 

In Fig. [14] we show the correlation between the different emis- 
sivities of our best solution at 30 GHz, where many emission 
mechanisms are present and the inversion is more challenging. 
These were obtained from the result of the jackknife tests pre- 
sented in Sect. l3.3l It can be seen that for almost all the templates 
the correlation is low. Three elements are to be highlighted. The 
first one is that the emissivities for the outer Galaxy molecular 
template, as well as the ionised template, both have pronounced 
departures from a Gaussian distribution. Looking at Fig.[T] both 
templates are patchy; the masking of several regions simultane- 
ously will change the morphology of each map significantly and 
therefore have a large impact on the derived emissivity. 

Secondly, rings 1-3 of the atomic and molecular phase are 
slightly anti-correlated (the contours are elliptical). This region 
of the Galaxy is indeed difficult to disentangle, and due to the 
low filling factor of these rings the algorithm has less to work 
with than in the solar circle, for example. Nevertheless the er- 
ror bars deduced from the jackknife analysis account for this 
slight degeneracy. Higher resolution maps will help reduce this 
crosstalk between templates in future analyses. 

Thirdly a small bias can be seen, as the results of the inver- 
sion are not always centred in the middle of the jackknife dis- 
tributions. We calculated this bias by looking at the difference 
between the mean deviation of each result with respect to the 
jackknife realisations, in units of result uncertainty. That is the 
bias B per template / and per band v can be written: 

B(i, y) = — i-lzf: [e(i, V) - ej(r, i, y)] , (10) 
cr(i, V) Ni- 

where A^, is the number of jackknife realisations, e(i,v) and 
cr(i, y) are the result and uncertainty of the inversion, and 
-ej(r,i,v) is the result of one jackknife realisation. We have 
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Fig. 14. Correlation matrix of the emissivities for each ring at 30 GHz. Overplotted in red are the results with corresponding error 
bars. One-dimensional histograms replace the contour plots where only one component is concerned. 



found this effect to be small, with a median bias of 0.005cr over 
all templates and frequencies, and with a mean absolute devia- 
tion of O.OSo". In all cases, the bias is always lower than 0.25cr. 



5.2. Effect of CMB subtraction 

The CMB removal in the Planck maps we have used may cre- 
ate significant residuals in the Galactic plane. Therefore, to 
test the impact of the CMB removal on our results, a different 
CMB estimate was reconstructed through a classical Internal 
Linear Combination (ILC) by means of Lagrange multipliers 
dEriksen et alj|2004l) . A ±15° strip in Galactic latitude was ex- 
tracted from the HFI CMB frequency channels maps (100, 143, 
217 and 353 GHz) reduced to a common resolution (10 arcmin), 
in units of Kcmb- This CMB component was subtracted from all 
the Planck LFl and HFI channels maps, resulting in a new and 
independent set of CMB-removed Planck maps. 

Even though the methods are similar, the resulting CMB 
in the Galactic plane is quite different. In constructing the the 
standard CMB subtrac tion (Needlet ILC, "NILC "), an intensity 
mask was applied (Pl anck HFI Core Teamll20lT ) which results 
in CMB residuals remaining in a substantial part of the Galactic 
plane. No such mask was used in the Galactic plane ILC. Outside 
of this mask, the differences between the two CMB estimates are 



of the same order as the residuals of our fit at low frequencies, 
44 GHz < V < 143 GHz. 

Using these maps to perform the same inversion, we have 
concluded that the impact is small, and almost always within the 
error bars. In Fig. [15] we plot the derived emissivities using the 
NILC and the Galactic plane centred CMB removal ("ILC"). 

To get a better idea of where the solution might be off by 
more than Icr, we show the difference in units of cr per band in 
Fig. [16] The value plotted is, for each component; 



6lLC ~ CNILC 



The agreement is generally very good, except for some 
points concentrated in the bands 143, 100 and 70 GHz. In these 
cases the CO templates are higher in the Galactic plane ILC 
whereas the H i and H^ templates are lower Nevertheless, there 
are no points that are different by more than 3cr. Therefore, 
these differences do not change our main conclusions, but would 
slightly modify the amount of anomalous emission, for example. 
As we are not modelling the specifics of any emission mecha- 
nisms at the lower frequencies, we are not sensitive to the impact 
of the CMB removal for this study. 
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Fig. 16. Comparison of emissivities derived after application of two different CMB subtractions. The value plotted is shown in 
Eq.[TT] The symbols represent Hi (plus), H2 (diamond), (triangle). Synchrotron (square) and dark gas (asterisk). Differences are 
generally small, but become larger in the range 143-44 GHz. 
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Fig. 15. Comparison of emissivities derived after application of 
two different CMB subtractions, the official Needlet ILC and the 
Galactic plane ILC. The atomic templates are the blue points, 
molecular are black, ionised is orange, dark gas is red and syn- 
chrotron is green. 
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Fig. 17. Dark gas intensity as a function of the residuals of the 
inversion at 30 GHz without (Left) and with (Right) the dark gas 
template. A distinct component of the residual (shown with the 
solid line) is seen to be correlated with the dark gas template. 
This correlated residual disappears when the dark gas is used in 
the inversion. 



5.3. Impact of using dark gas template 

The dark gas template used in this study has been calculated by 
iPlanck Collaboration! (|201 id) . It is the result of finding at what 
column densities the correlation between total hydrogen column 
density (Hi+2H 2) versus optical depth of the dust departs from 
linearity. The frequencies used to obtain the map at 857 GHz that 
we use as a template are 857 and 545 GHz, along with the IRAS 
100 i-im band. There is some circularity in using a template based 
on a residual in the Planck data to fit the same data. However, the 
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data used is restricted to the two highest HFI frequencies. Here 
we present the impact at 30 GHz, which in no respect suffers 
from the same circularity, as the data are far from the frequencies 
used to derive the dark gas template, and they come from LFI, a 
completely separate instrument from HFI. 

We have performed the inversion with and without the dark 
gas template at 30 GHz, and looked at the correlation between 
the residuals and the column density of the dark gas. The re- 
sults are shown in Fig. [TtI the inversion performed without the 
dark gas template is on the left, the same inversion with the dark 
gas template is on the right. A distinct component of the residu- 
als can be seen to correlate very well with the dark gas column 
density when the dark gas template is not used. This correlation 
completely disappears when it is used. In addition, the reduced 
for the fit shows a significant improvement of approximately 
50%. 

The choice of both the opacity correction for the H i and the 
XcQ value chosen will have an impact on our results. We have 
chosen to keep the analysis simple for this study, as the spin tem- 
perature we assume, together with the dark gas template, should 
account for most of the optically thick gas in the plane. The abil- 
ity of the dark gas template used to adequately account for all the 
optically thick gas remains to be studied, so future work could 
involve using more sophisticated methods to account for this ef- 
fect. 



6. Conclusions 

The analysis performed in this study provides a reaUstic descrip- 
tion of the dust and gas properties as a function of Galactic ra- 
dius in the Milky Way. The dust temperature in the H i is seen 
to decrease as a function of Galactocentric distance from 24.0 to 
13.9 K, and the temperature in the molecular phase is heated by 
star formation associated to the Galactic spiral arms. The opac- 
ity found for grains in the solar circle is t/A^h - 0.92 + 0.05 x 
10"^^ cm^, and with no significant variation with Galactic radius, 
even though the dust temperature is seen to drop by over 10 K 
from the centre to the outer Galaxy. The dust temperature in the 
Hi gas in the solar circle, 17.6 + 0.1 K is also compatible with 
the recent Planck value of high-latitude cirrus (17.8 ± 0.9 K). 

The extension of our analysis to lower frequencies than pre- 
viously used in inversion techniques has allowed us to place con- 
straints on the free-free, synchrotron and anomalous microwave 
emission: 

- Anomalous dust emission is clearly seen in the atomic, 
molecular and dark gas phases. It is well fit by a very simple 
model consisting of spinning PAH molecules embedded in 
each of the gas phases. We highlight regions where spinning 
dust emission may be strong with respect to other emission 
mechanisms. According to our simple model, in the Galactic 
plane spinning dust accounts for 25 + 5%(statistical) of the 
total emission at 30 GHz. However, systematic uncertainties 
linked to our model may make this value more uncertain. 

- The dark gas phase has been explored spectrally across the 
Galactic plane for the first time. Its SED is similar to the 
molecular phase and therefore this phase seems to be tracing 
diffuse H2 where CO is disassociated by energetic photons, 
as well as tracing optically thick CO and H i. 

- The free-free emission is completely traced by our ionised 
component and indicates that the free-free estimate from 
WMAP component separation methods may be too high. In 
the ionised phase, spinning dust may be present but it is dom- 
inated by free-free emission. 



- Synchrotron emission is well characterised by a power law 
with a spectral index of - 1 .0. 

This first analysis of the large-scale Galactic emission as 
seen by Planck has constrained the large scale ISM properties in 
the Galaxy. It could be improved in the future via several modi- 
fications. Individual regions could be studied using higher angu- 
lar resolution data, which would probe smaller scale fluctuations 
that have been missed by the present study. The radial separa- 
tion could be done differently to explore the difference between 
arm and inter-arm ISM, in order to focus more on the role of 
ISM conditions on the star formation process. It should also be 
possible to obtain an independent estimate of the dark gas in the 
plane, and ideally some radial information about its behaviour. 
In the near future, Planck component separation products could 
be used in the place of WMAP ones. Use of radio recombination 
line surveys currently underway would be very useful to con- 
strain the free-free emission in the plane, and as a function of 
Galactic radius. 
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Table A.l. Result of the inversion per band and per component in MJy sr~^ / 10^" H atoms, except for the synchrotron component 
in which case it is in MJy sr"^(v) per MJy sr"^(408 MHz). 



Band [GHz] 


Hi(l) 


Hl(2) 


Hi(3) 


Hl(4) 


Hi{5) 


25000 


1.4 ±0.4 xlO"' 


5 ±3 xlO-2 


7.8 ±0.6x10"^ 


2.4 ±0.2x10-2 


6.0 ±0.9x10"' 


12000 


<1.5 xlO"' 


<4.6 xlO-2 


8.2 ±1.3 xlO"2 


<1.2 xlO"2 


<5.9 XlO"' 


5000 


8 ±3 xlO"' 


3 ±2 XlO"' 


4.1 ±0.5 XlO"' 


<6.2 xlO"2 


<2.5 xl0"2 


3000 


3.3 ±1.5 


2 ±1 


1.8 ±0.2 


5.9 ±0.6x10"' 


<4.6 xl0"2 


2143 


4.9 ±1.2 


3.0 ±0.8 


2.6 ±0.2 


1.3 ±0.1 


2.6 ±0.3 XlO-' 


1250 


1.5 ±0.6 


1.6 ±0.4 


1.5 ±0.1 


9.7 ±0.4 XlO-' 


4.0 ±0.1 XlO-' 


857 


8 ±2 xlO-' 


9.2±1.7xl0-' 


7.7 ±0.3 XlO-' 


5.7 ±0.1 XlO-' 


3.2 ±0.1 XlO-' 


545 


1.9 ±0.7 xlO-' 


2.9 ±0.5 XlO"' 


2.3 ±0.1 XlO-' 


2.0 ±0.1 XlO-' 


1.4 ±0.1 XlO-' 


353 


3.8 ±1.5 xlO-2 


6.7 ±1.2x10^2 


5.4 ±0.2x10-2 


4.9 ±0.1 xlO-2 


4.0 ±0.1 xlO-2 


217 


6 ±3 xlQ-^ 


1.4 ±0.3 xl0"2 


1.1 ±0.1 xl0"2 


1.1 ±0.1 xlO-2 


9.5 ±0.1 xlO"' 


143 


1.1 ±0.7 xlO-2 


2.7 ±0.5 XlO"' 


2.4 ±0.1 XlO"' 


2.3 ±0.1 XlO"' 


2.3 ±0.1 XlO"' 


100 


<7.3 xlO"" 


6 ±2 xlO"* 


5.8 ±0.6 XlO"'' 


5.9 ±0.3 XlO"" 


7.6 ±0.1 XlO"" 


70 


<4.1 xlO"" 


1.1 ±0.9 xlO"'* 


1.6 ±0.5 XlO"" 


1.3 ±0.1 xlO"" 


2.9 ±0.1 xlO"" 


44 


<3.3 xlO"" 


<8.9 xlO"' 


9 ±4 xlQ-' 


8.3 ±1.3 xlO-' 


1.4 ±0.2x10-" 


30 


<3.4 xlO " 


1.5 ±1.2x10"'* 


1.3 ±0.4 X 10-* 


1.4 ±0.1 xlO-* 


1.2 ±0.2x10-" 


1.4 


<1.1 xlO"^ 


<3.8 xlO"'* 


<3.1 xlO"* 


<1.7 xlO-* 


4 ±2 XlO"" 


Band [GHz] 


H,(l) 


H,(2) 


H,(3) 


H,(4) 


H,(5) 


25000 


2.3 ±1.3 xlO"^ 


6.2 ±0.6x10"^ 


9 ±6 XlO"' 


8 ±4 XlO"' 


3.0 ±0.6xl0"2 


12000 


4.8 ±1.9 xlO-2 


5.1 ±1.0x10-- 


4.8 ±1.1 xlO"2 


<1.8 xlO-2 


<3.6 xlO-2 


5000 


1.7 ±1.3 xlO"' 


5.7 ±0.7x10-' 


2.5 ±0.7x10"' 


<8.1 xlO-2 


<1.8 XlO-' 


3000 


1.4 ±0.6 


2.3 ±0.3 


8 ±2 XlO"' 


2.7 ±1.0x10-' 


9 ±2 XlO-' 


2143 


2.8 ±0.4 


3.6 ±0.2 


1.8 ±0.2 


4.3 ±0.9 X 10-' 


1.5 ±0.3 xlO" 


1250 


2.0 ±0.2 


2.1 ±0.1 


1.3 ±0.1 


5.3 ±0.4x10-' 


1.1 ±0.1 XlO" 


857 


8.0 ±0.8 xlO"' 


8.3 ±0.4x10"' 


6.1 ±0.4x10"' 


3.8 ±0.2 XlO"' 


6.1 ±0.5 XlO"' 


545 


2.1 ±0.2 xlO-i 


2.1 ±0.1 XlO"' 


1.7 ±0.1 XlO"' 


1.4 ±0.1 XlO"' 


2.1 ±0.1 XlO"' 


353 


4.3 ±0.5 xlO-2 


4.1 ±0.3x10"^ 


3.7 ±0.2xl0"2 


3.7 ±0.1 xlO"2 


5.3 ±0.3 xl0"2 


217 


1.1 ±0.1 xlO"^ 


1.0 ±0.1 xlO"2 


9.5 ±0.4 xlO"' 


1.0 ±0.1 xlO"2 


1.4 ±0.1 xl0"2 


143 


1.4 ±0.2 X 10"^ 


1.3 ±0.1 xlO-' 


1.2 ±0.1 xlO-' 


1.7 ±0.1 xlO-' 


2.3 ±0.2 X 10-' 


100 


1.9 ±0.1 xlO"^ 


1.6 ±0.1 xlO-' 


1.5 ±0.1 xlO-' 


1.6 ±0.1 xlO-' 


1.7 ±0.2 xlO-' 


70 


2.1 ±0.5 xlO"" 


1.8 ±0.2 xlO-'' 


2.2 ±0.3 xlO-* 


2.3 ±0.1 xlO-* 


3.7 ±1.6 X 10-" 


44 


1.9 ±0.4 xlO"'' 


1.7 ±0.3 xlO-'' 


2.4 ±0.3 xlO-* 


1.3 ±0.1 xlO"* 


4.4 ±1.9 xlO-" 


30 


2.8 ±0.4 xlO"" 


2.5 ±0.4 xlO"'* 


3.1 ±0.3 xlO"* 


1.3 ±0.1 xlO"" 


6 ±2 xlO"" 


1.4 


<3.6 xlO"" 


<1.6 xlO"'* 


<1.7 xlO"* 


<8.5 xlO"' 


<5.9 xlO"' 


Band [GHz] . 


Hii 


Sync. 


Dark 






25000 


9.7 ±0.9 xlO"' 




4 ±3 XlO-' 






12000 


3.3 ±0.3 




<1.5 xlO-2 






5000 


1.9 ±0.1 xlO' 




<7.1 xlO"2 






3000 


2.8 ±0.1 xlO' 




<1.7 XlO"' 






2143 


2.9 ±0.1 XlO' 




1.0 ±0.1 






1250 


1.2 ±0.1 XlO' 




8.3 ±0.4x10-' 






857 


4.7 ±0.2 


- 


4.4 ±0.2 XlO-' 






545 


1.3 ±0.1 




1.5 ±0.1 xlO ' 






353 


3.1 ±0.1 xlO^' 




3.5 ±0.1 xlO-2 






217 


7.9 ±0.2x10-2 




7.9 ±0.2 XlO-' 






143 


3.0 ±0.1 xlQ-^ 




1.8 ±0.1 XlO-' 






100 


2.3 ±0.1 xlO-2 


3 ±2 xlO"3 


7.1 ±0.3x10-* 






70 


2.2 ±0.1 xlO-2 


<4.7 XlO"' 


2.5 ±0.2x10"* 






44 


2.3 ±0.1 xlO-2 


9 ±3 XlO"' 


1.5 ±0.2 XlO"'' 






30 


2.2 ±0.1 xlO-2 


2.0 ±0.3 xlO"2 


1.7 ±0.2x10"" 






1.4 


1.8 ±0.1 xlO-2 


2.9 ±0.1 XlO"' 


<3.5 xl0"5 
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